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Abstract 

In the framework of renormalization-group improved cosmologies, we analyze both theoretically 
and observationally the exact and general solution of the matter-dominated cosmological equations, 
using the expression of A = A(G) already determined by the integration method employed in a 
previous paper. A rough comparison between such a model and the concordance ACDM model 
as to the magnitude-redshift relationship has been already done, without showing any appreciable 
differences. We here perform a more refined study of how astrophysical data (Union2 set) on type-I 
supernovae, gamma ray bursts (in a sample calibrated in a model independent way with the Snela 
dataset), and gas fraction in galaxy clusters (using a sample of Chandra measurements of the X-ray 
gas mass fraction) affect the model and constrain its parameters. We also apply a cosmographic 
approach to our cosmological model and estimate the cosmographic parameters by fitting both the 
supernovae and the gamma ray bursts datasets. We show that this matter-dominated cosmological 
model with variable Newton parameter and variable cosmological term is indeed compatible with 
the observations above (on type la supernovae, the gamma ray bursts Hubble diagram, and the gas 
mass fraction in X-ray luminous galaxy clusters). The cosmographic approach adopted confirms 
such conclusions. Finally, it seems possible to include radiation into the model, since numerical 
integration of the equations derived by the presence of both radiation and matter shows that, after 
inflation, the total density parameter is initially dominated by the radiation contribution and later 
by the matter one. 
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INTRODUCTION 



Pushed on by the overwhelming flow of observational data in the last fifteen years, most 
cosmologists today agree on a well defined cosmological paradigm, based on General Rela- 
tivity plus a cosmological constant A. This paradigm is known as the Concordance Cosmo- 
logical Model and accounts not only for the early formation of large-scale structures but 
also for the more recently discovered stage of acceleration of the universe. As to the building 
up of galaxies and galaxy clusters, it has been necessary to introduce an ingredient like dark 
matter, which was first employed to succeed in describing rotational curves in spiral galaxies 
[2! (for an alternative view, see for example the work in Ref. jsj). On the other hand, cosmic 
acceleration requires the consideration of the so-called dark energy as the major ingredient 
of the cosmic content ^|, A being just the simplest way to consider it. Both of them sum 
up to more than 95% of the matter-energy constituents around us. 

Such dark energies, on the other hand, have simply hidden the fundamental issues, since 
so far no exhaustive physical explanation for them has been put on firm theoretical and 
experimental ground. This has led to many alternative ways to reproduce the astrophysical 
phenomena cited above, not only by introducing theoretically well motivated new particles 
and fields, but also by means, for instance, of possible geometrical changes of the spacetime 
structure (/(-R) theories are well known examples of this kind of proposals jsl 

Without entering any details of so many attempts to describe the cosmological behaviour, 
we shall here consider only some aspects of one of them, i.e. that stemming from the 
possibility that not only the cosmological term A could vary with space and time, but also 
the gravitational coupling G. In this context, one way to achieve the physical realization of 
such assumptions is to study cosmological dynamics by analyzing "renormalization group 
(RG) induced" quantum effects, which drive the dimensionless cosmological "constant" Mk ) 



and Newton "constant" g{k) from an ultraviolet attractive fixed point 



Bis 



9( M [nl M- 



In the exact theory, such a non-Gaussian ultraviolet fixed point implies its nonperturbative 



renormalizabihty [l2|[l3|Q[l5|[l6|[l7|[l8|[l9||20^ result, this Renormalization 

Group -improved framework describes gravity at a typical distance scale i =Ji~^, and 
introduces an effective average action rfc[(7^;^] for Euclidean quantum gravity [8], finally 
implying an exact functional Renormalization Group equation for the fc-dependence of F^. 



This framework is usually known as quantum Einstein gravity. Within it, one can get 
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an explicit /c-dependence of both the running Newton and cosmological terms G{k) and 



A(A;), which can 
smg 



)e relevant both for the initial Planck era and the structure of black hole 
23 



23|. 



Taking into account its inherent infrared divergences, quantum Einstein gravity can be 
subject to strong renormalization effects even at very large distances. In cosmology, such 



effects lead to a dynamica 
logical constant problem 



relaxation of A and can also be assumed to deal with the cosmo- 



24l |. Viewing the late accelerated expansion of the universe as a 



renormalization group evolution near a non-Gaussian infrared fixed point 25| (although the 
actual existence of an infrared fixed point has not been proved as yet), one can assume that 
the transition between standard FLRW cosmology and accelerated Renormalization Group 
driven expansion occurs at the time when the fixed point is almost reached. As to this, some 
agreement has been found between this kind of model and early SNela observations 26|. 



As a matter of fact, for a homogeneous and iso trop ic universe, it is possible to identify k 
with the inverse of cosmological time, k (xl/t 2l|] 25|], hence deriving a dynamical evolution 
for G{k) and A(fc) induced b y th eir Renormalization Group running. The Arnowitt-Deser- 
Misner (ADM) formulation 271] builds a modified action functional which reduces to the 
Einstein-Hilbert action when G is constant. Within such a framework, and always assuming 
homogeneity and isotropy, one can obtain a power-law growth of the scale factor for pure 
gravity and for a massless theory, in agreement with what is known on fixed-point 
cosmology. On the other hand, by means of the so-called Noether Symmetry Approach 28] 



29(1, in Ref. 30|| we have also proposed solutions for the pure gravity case, which mimic 



inflation without introducing a scalar field in the cosmic content. In Ref. |31|, for gravity 
with a scalar field, this approach only succeeds in fixing the expressions of A(G) and V{ip) as 
A oc G and V cc ip^, respectively, while the transformed cosmological equations derived by 
means of the method do not seem to be easily solvable 3l|], therefore giving no new insight 
into possible solutions. 

In what follows, we take again into account the exact solutions of the fiat dust matter- 
dominated cosmological equations (without any scalar field), already investigated in Ref. [s^ 
by means of the Noether Symmetry Approach, and by using an expression for A = A(G') 
determined by the method itself. After briefly reviewing the theoretical model, we show that 
our cosmological model is compatible with various recent observational data, in particular 
with the observations of type la supernovae (SNela) (we use the recently updated SNela 
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sample, referred to as Union2 33|], containing 557 SNela spanning the redshift range 0.015 < 
z < 1.55.), the Gamma Ray Bursts Hubble diagram (GRBs HD) (we use a sample calibrated 
in a model independent way with the Snela dataset [sJ]), and the gas mass fraction in X- 
ray luminous galaxy clusters (we use a sample of Chandra measurements of the X-ray gas 
mass fraction in 42 hot {kT > 5keV), X-ray luminous, dynamically relaxed galaxy clusters 
spanning the redshift range 0.05 < z < 1.1 35l|). 

We are aware that such observational tests cannot definitively distinguish our model from 
the standard ACDM one, at least because of the presence of still large observational errors. 
However, it provides a crucial and necessary test of reliability. 

Looking at future studies on this, we finally apply to our cosmological model a cosmo- 
graphic approach, which can indeed contribute to select realistic models without imposing 
arbitrary choices a priori. As a matter of fact, cosmography and its reliability are based 
on the assumptions that the universe is spatially homogeneous and isotropic on large scale, 
and luminosity distance can be "tracked" by the derivative series of the scale factor a{t). 
We actually estimate the cosmographic parameters here derived by fitting both the SNela 
Union2 dataset and the calibrated GRBs HD. 

We also begin to study how our model is affected by the inclusion of radiation into 
the cosmic content. It in fact turns out that, by performing the numerical integration 
of the equations so rewritten, the total density parameter is initially dominated by the 
radiation contribution and later on by the matter one, leaving then space to the now observed 
accelerated stage. 

The scheme of the paper is as follows. In Section 2 we summarize the Lagrangian formu- 
lation used to derive the Renormalization-Group improved Einstein cosmological equations 
with the ordinary matter energy-momentum tensor, as well as the results deduced from the 
Noether symmetry found in Ref. j32|. In Sections 3, 4 and 5 we present the comparison 
of theoretical predictions with observational data. Section 6 is then devoted to the cos- 
mographic approach, and Section 7 to inclusion of radiation into the model. Finally, some 
conclusions are drawn in Section 8. 
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THEORETICAL MODEL 



Let us consider the approach outhned in Ref. [27| [30| [32| [3l| and there apphed to models 
of gravity with variable G and A in the context of quantum Einstein gravity. It is known 
that, in a homogeneous and isotropic universe, an independent dynamical G is equivalent to 
metric-scalar gravity already at classical level 36| 37|], while independent variations (with 



position and time) of G and A can lead to pathological situations. Indeed, if A were an 
independent variable, one should write that the momentum conjugate to it vanishes, and 
the preservation in time of this primary constraint would imply a vanishing lapse function 
and hence a "collapse" of spacetir ne g eometry |271] . All this in fact leads to assume a generic 
functional dependence A = A(G) |27l |. 

In the matter-dominated case in a flat homogeneous and isotropic cosmology (with a 
signature —,+,+,+ for the metric, lapse function = 1 and shift vector A^* = 0), as in 



Ref. 



32| , we start from the Lagrangian 
1 



1 



SttG 



3aa — a A H — ^a: 



.G^ 



2' G2 



(1) 



where G = G{t) and A = A{G{t)), while dots indicate time derivatives, and u i s a nonva- 
nishing interaction parameter introduced in Ref. {27 1 and also used in Ref. 32 1, where it 
was shown that fi > 2. The matter contribution is of course given by = —Da~^^"'~^\ 
with 1 < 7 < 2; here, we have to take 7 = 1 for dust, while D is a suitable integration 



constant connected to the matter content. As shown in Ref. 
Euler-Lagrange equations for a and G 



32|, from Eq. ([T]) we get the 



• 2 

a a 
a 20? 



A 
2" 



aG fiG^ 
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3 G^ 



G 



The Hamiltonian constraint 



32| 



d2 



fiG^ 
QG^ 



StcG 



Da 







(2) 
(3) 

(4) 



a" a Lr" 3 

is equivalent to the constraint on the energy function associated with L [28|[29| [30| |32| [31 1 



For matter, the dust case involves a zero pressure, pm = 0, and an energy density pm = 
Da^^, so that the matter term in the Lagrangian is simply a constant. It therefore has no 
effect on the equations of motion with respect to the pure gravity case; nevertheless, it has 
to be considered, since it occurs in the constraint equation The system of equations 
of motion can then be solved ^] by using the Noether Symmetry Approach 28| 29|], in 
which we consider L as a point Lagrangian, a function of the variables a and G, and their 



first derivatives 



29| 30| |32|. We have already shown that a consistent choice of the 



function A = A(G') leads to the existence of a Noether symmetry for the Lagrangian 32 |. 
As a matter of fact, in the matter-dominated case we get the same Noether symmetry as 
in the pure gravity situation, so that, by using the same transformations introduced in this 
latter case, one can write a = a{t) and G = G{t) as therein, now just updating the energy 
constraint. (For more details, see Ref. (32].) 

The dynamics of A is coupled with that of G and is driven by the equation 

(6) 



2(1-.)A + G§ 



0, 



where the parameter J is an arbitrary constant related to the interaction factor p by the 
relation fi = |(3 — 2J)^ 7^ 0, |. This equation admits the solution 

A = A{t; n) = WG^ , (7) 

W being an integration constant and n{J) = ^x^jy- K turns out that W is related to the 
present value of G |32|, and we might determine W in order to get Gq = Gj\f = GT^ton- 



Furthermore, we fix time scale and origin so as to get a(0) = 0. Thus, we find |32| 



a = a(t) 
G = G{t) 
where we define the constants 
A 



A t6 



B + t6n-r 



3n-l 



(8) 
(9) 



n(l + 6n) 12-n? 12n'^ 

A{n, W) = 12^=^n^{Qn - l)^{12n - l^W 



B = B{n,W,D) = W~^ 2'-^^=^{3n)^{6n-l)'^{12n-l)7iD , 
C = C{n, W) = {6n - 1)2(3— i) [12^2(12^ - 1)] IV^n-i _ 

Here we point out that the asymptotic time behavior of the scale factor is characterized by 
the two exponents 



(10) 

(11) 
(12) 



Pi = 



12n^ 
(6n - 1) 



and we want that pi > 1 and p2 < 1 (in order to obtain early matter domination and a 
later accelerated evolution ), which implies that we find a limited range of variability for the 
n parameter jsol], (3 — -\/3)/6 < n < {3 + y/3)/Q. We then see that, when D = 0, one has 
B = 0, and we recover the same results obtained in Ref. [sij for the pure gravity model. 
Eqs. fIlOp and fill I) make it possible to obtain D and W as functions of A, B and n, i.e. 

2 Q 1 1 I -I 1 1 

_ 26"-i 36"-in6"-i^ (6n — 1) 1-6" 

= — , \^^} 

71 

W = 12^+V^(6n- l)^(12n- 1)A^ . (15) 

To make things analytically simpler and obtain a better control of the space of parameters 
related to the integration constants, we set the present time to = 1. This fixes the scale of 
time according to the (unknown) age of the universe. In other words, this means that we 
are using the age of the universe, to, as a unit of time, and the whole history of the universe 
has been squeezed to the range of time [0, 1]. We then set = cf(l) = Ij which is standard, 
and finally Hq = H{1) 1. Because of our choice of time unit, it turns out that our Hq is 
not the same as the Hubble constant which appears in the standard FLRW model. 
Such choices introduce a constraint between A and B: 

A={l + By. (16) 

On choosing to normalize a^, the definition of the redshift z^oq/o — l = l/a— 1 yields 

z = z{t) = {B + lyt^^ [b + 1 6;^+^) " - 1 , (17) 

depending only on the parameter n. It turns out that we obtain the following expression for 
Hq and Aq: 

i/„^J/fa^l)=, «P + ^>"^ (18) 

A„.Afa^l)^- °(« + y^+l)"-' . (19) 

on — 1 

We have constrained the parameters of our model in order to have Go = Gn = 1- Therefore, 
the density parameter of matter is 



^ 87rG{t,B,n)D ^ SnGoD _ B{B + l){6n-l) 
™" 3Hit,B,nf 3Ho' ~ %B + 2Yn^ ' ^^^> 
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and it turns out that the constraint 



rim + ^AO + ^GO = 1 (21) 

is satisfied. Moreover, as aheady pointed out in Ref. [38|, we do not expect that Gq = Gn = 
(jNewton, evcu if we fixed the parameters in order to get it. Anyway, we have to consider that 
small differences between Go and Gn could imply Go varying a lot in time, with relevant 
effects on the evolution of the universe. It is worth noting that G has a special role in 
the subject of time variation of the fundamental parameters. Actually a dependence of G 
on time may point out to violation of the strong equivalence principle, but not necessarily 
of the Einstein equivalence principle, whereas the nonconstancy of the other "constants", 
like the electroweak or strong coupling constants, necessarily represents a violation of the 
equivalence principle in both its forms. Here we will only discuss some aspects of the response 
of primordial abundances to the time variations of G, as expected in our model, on the basis 



39| 



of previous analysis on the subject performed in Ref. 

It turns out that the production of each element responds in its own way to a variation 
6G of the Newton constant. A general study, that can account for a time dependence of G 
during the BBN period, requires the introduction of suitable functions which describe the 
response of each element to an arbitrary time-dependent modification of the early universe 
expansion rate, and it is out of the aim of the present work. So, we simply discuss the 
observational bounds on the possible variations of the gravitational constant in the early 
universe, considering the best limit (at 3cr), 6G = 0.09^q^I, obtained in Ref. Sol, by 
combining observational results with the measurements of the baryon to photon ratio 
obtained from CMB and LSS data 5^]. We will postpone to a forthcoming paper the detailed 
analysis of the dependence of the various elemental abundances on the time variation of G 
for our model. It turns out that our model can satisfy such a best fit limit, provided that 
the n parameter is appropriately selected (the role of B is only marginal with respect to 
this strong constraint), and compatible with the other basic cosmological observations, as 
shown in Fig. [H 

To further investigate this issue, we first evaluate the fractional time rate of change of G 

^~xlO-^ (22) 
Go 

Here, it is important to remember that we are using the age of the universe as unit, so that 



0.004 - 



-0.004 - 



-0.006 - 




FIG. 1: Time evolution for the relative variation 6G for our model, with B = 2.77, n = 0.32. As 
we will discuss in the next sections such values are fully consistent with the constraints resulting 
from other cosmological datasets. 

the effective rate is of order 10~^^ yr^^, in agreement with the observations (see for instance 



Ref. 



We can thus say that, eventually, the running of both the gravitational coupling G = G{t) 
and the cosmological term A = A(t) induced by quantum effects appears to yield both a 
primordial inflation (soon after the universe exits the region where the attraction basin of 
the non-Gaussian (ultraviolet) fixed point works [38j]), and a later inflationary epoch in a 
matter- dominated period of the expansion of the universe. This appears to be interesting 
since it is obtained always without having to introduce any scalar field in the cosmic content 



32|. 



Re-parametrization of the model 

Let us now exhibit a re-parametrization of our model in terms of Hq and Vim instead of 
B and n. On inverting the systems of Eqs. f|T8l) and fl20|) it is indeed possible to recover the 
parameters B and n as functions of Hq and VLm'-, actually, we have that 

2 — 3HoQm 



1 ( 

n = ■ 

24 V 



(^VS^ HoiSHorim + 2) (9iJ2fim + GHo - 8) + 9H^n^ + 6Ho^ . (24) 
The conditions on the two exponents pi and p2 (that is pi > 1 and p2 < I) give rise to a 
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FIG. 2: Allowed regions (in light blue) of the space of parameters n and B. 



constraint on the space of parameters for Hq and Qm, as shown in Fig. |2l actually, it turns 
out that 

2 < i/o(3i^o^^m + 2) < 4. (25) 

It is worth noting that the space of parameters Hq and Qm is reasonably deducible from 
physical arguments. As a matter of fact, Q^n varies in the range [0, 1], and the range of 
variation for Hq can be inferred by assuming that the age of the universe can be written in 
the following way: 

to = 7 X IGy = 3.15 ■ 10^''7s, (26) 

where 7 is a constant to be determined by astronomical observations. With this definition 
the value of Hq can be related to the small h = -ffo/lOO of the standard FLRW model. It 
turns out that 

Ho = 0.1 h-f. (27) 



If we accept that to = 13.76 ± 0.11 and h = 0.71 ±0.014, as given by WMAP7 |4l|, then 
the region of variability at 2a for Hq turns out to be (0.92 , 1.03). It is interesting to note that 
the constraint in Eq. (l25l) is satisfied in the whole domain of and for 0.06 < Qm < 0.6. 



CONSTRAINTS FROM RECENT SNEIA OBSERVATIONS 

Over the last years the confidence in type la supernovae as standard candles has been 
steadily growing. Actually, it was just the SNela observations that gave the first strong 
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indication of an accelerating expansion of the universe, which can be accounted fo r by 



assuming the existence of some kind of dark energy or nonzero cosmological constant 



42|. 



Since 1995 two teams of astronomers - the High-Z Supernova Search Team and the Supernova 



Cosmology Project - have been discovering type 



of both teams were published in Refs. 



421 and 



a supernovae at high redshifts. First results 



43|. 



As to a first comparison od theory with observations, we here consider the recently 



updated Supernovae Cosmology Project Union2 compilation [33|, which is an update of 
the original Union compilation, now bringing together data for 719 SNela, drawn from 
17 datasets. Of these, 557 SNela, spanning the redshift range (0.015 < z < 1.55.), pass 
usability cuts and outliers removal, and form the final sample used to constrain our model. 
We actually compare the theoretically predicted distance modulus /i(z) with the observed 
one, through a Bayesian approach, based on the use, as merit function, of the likelihood 
C = exp (— |x^)- The distance modulus is defined by 

^ = m-M = 51ogc/i(2;) + 51og(^) +25, (28) 

where m is the appropriately corrected apparent magnitude including reddening, K correc- 
tion etc., M is the corresponding absolute magnitude, and d^ is the luminosity distance in 
Mpc. However, in our cosmological model with variable G and A, it is important also to 
include in Eq. ( |28|) corrections describing the effect of the time variation of the gravitational 
constant G on the luminosity of high redshift supernovae. If the local value of G at the 
spacetime position of the most distant supernovae differs from Gn, this could in principle 
induce a change in the Chandrasekhar mass Mch oc G~^. Some analytical models of the 
supernovae light curves predict that the peak luminosity is proportional to the mass of nickel 
produced during the explosion, which is a fraction of the Chandrasekhar mass. The actual 
fraction varies in different scenarios, but the physical mechanism of type la supernovae ex- 
plosion always relates the energy yield to the Chadrasekhar mass. Assuming that the same 
mechanism for the ignition and the propagation of the burning front is valid for SNela at 
high and low redshifts, the predicted apparent magnitude will be fainter by a quantity [44 1 

AMa = ^ log (^] . (29) 



4 "\Go, 

Taking this into account the distance modulus becomes 



m- M = 5\ogdL{z) +5 log (j^) + 25 + AMg . (30) 
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The presence of this correction allows us to appropriately test our model by using the SNela 



sample 



4| |45(|46|. 



In our flat and homogeneous cosmological model the luminosity distance can be expressed 
as an integral of the Hubble function as follows: 

1 



(31) 



where H{z) is the Hubble function expressed in terms of redshift. It turns out that the 
luminosity distance can be expressed as a function of time in the following way: 



6{B + 2)n'^{t -1){B + (t^ [b + 



Qn — 1 



1 1 
-n + 1 - — , n; -n + 2 - — ; 

bn bn 



6n \ n 
B + t 6n-l \ / 
X I — 1 (t6"-l (i?+t6n-l 



6n 

1 1 t6"-l 

-n+l- —,n]-n + 2- —] — 

bn bn B 




where 2-^1 is an hypergeometric function. On inverting the relation z{t) in Eq. ( ITTI) . we 
obtain 



= 2^-W((z + l)(fi + l)-") 



-l/n 



(33) 



xW((^ + l)(5 + l)-)- LB2((^+ 1)-")- +4 -B 



and we can construct di^z) and evaluate the distance modulus according to Eq. fl30|) . so as 
to perform our likelihood analysis, by maximizing the likelihood C = exp (— on a grid 
in the space of parameters 58j B and n. In order to constrain the parameters of our model 



only, when we perform our statistical analysis with the Union2 compilation, we marginalize 
over h, that is, we maximize the likelihood Cmarg = j^""^^ dhexp (— where hmin and 
hmax are fixed by using the latest WMAP7 results. We obtain xleduced ~ 0-97 ^^i 557 data- 
points and the regions of confidence at 3 a for Hq and Qm are (0.92,1.01) and (0.24,0.4), 
respectively. If we do not marginalize over h, we obtain hhest = 0.70^*^0 02, from which we 
can infer the following interval of confidence at 3(t: r G (12.8, 14.7) Gyr. In Fig. [3] we plot 
the Union2 data set with the best fit modulus of distance, showing that they are, indeed, 
well-fitted by our model. 
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FIG. 3: The fit of the Union2 data set with the theoretical modulus of distance fi with respect to 
redshift z. 

CONSTRAINTS FROM CALIBRATED GAMMA RAY BURSTS HUBBLE DIA- 
GRAM 

Let us now assess another set of observations which we think essential to begin to un- 
derstand how much our theoretical model can be considered as a reliable one. As a matter 
of fact, it has been recently empirically established that some of the directly observed pa- 
rameters of Gamma Ray Bursts are connected with the isotropic absolute luminosity Liso, 
the collimation corrected energy Ej, or the isotropic bolometric energy Eiso of a GRB. Such 
observable properties of the GRBs include the peak energy, denoted by E^ i, which is the 
photon energy at which the u F^, spectrum is brightest; the jet opening angle, denoted by 
^jet, which is the rest-frame time of the achromatic break in the light curve of an afterglow; 
the time lag, denoted by Tiag, which measures the time offset between high and low energy 
GRB photons arriving on Earth; and the variability, denoted by V, which is the measure- 
ment of the spikiness or smoothness of the GRB light curve. In the literature, there is a 
wide variety of choices for the definition of V ^] in which the observed V value varies as 
the inverse of the time stretching, so the corresponding measured value should be multiplied 
)y a correcting factor {1 + z). An additional luminosity indicator is the minimum rise time 



471] denoted by Trt, and taken to be the shortest time over which the light curve rises by 
half the peak flux of the pulse. 

These quantities appear to correlate with the GRB isotropic luminosity, its total 
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collimation-corrected or its isotropic energy. This property cannot be measured directly 
but rather it can be obtained through the knowledge of either the bolo-metric peak flux, 
denoted by Pboio, or the bolo-metric fluence, denoted by Shoio- Therefore, the isotropic 
luminosity is given by 

Liso = 47r4(z)Pboio , (34) 

the total isotropic energy reads as 

Eiso = 47r4(z)5boio(l + z)-' , (35) 
and the total collimation-corrected energy is 

= 47r4(z)5boloi^bcam(l + (36) 

where -Fbeam is the beaming factor. The correlation relations are power-law relations of either 
Liso or E^ or E^^o as a function of riag, V, -Epeak, ^rt, i.e. 

Tp U r^'^iso.pcak 

-C/iso — Oiso,peak-C/pcak ' 



77" t. r^"7,peafc 

— (^7,peafc-C/peak ' 

r I, rpO-peak 

— OpeakJ^peak ■ 

(37) 

Therefore, Lj^o, E^ and i^iso depend not only on the GRB observables Pboio or Sboio, but also 
on the cosmological parameters, through the luminosity distance d^lz). As a consequence, 
there is a flerce problem to overcome, since it is not immediately possible to calibrate such 
GRBs empirical laws, and to build up a new GRBs Hubble diagram, without assuming any 
a priori cosmological model (which is known as the circularity problem). 

In Refs. {mi and {48 1 we have applied a local regression technique to estimate, in a model 
independent way, the distance modulus from the recently updated Union SNela sample, con- 
taining 557 SNela spanning the redshift range of 0.015 < z < 1.55. The derived calibration 
parameters have been used to construct an updated GRBs Hubble diagram. In particular, 
by using such a technique, we have fltted the so-called Amati relation (see Ref. [34] for 
details) and constructed an updated Gamma Ray Bursts Hubble diagram, which we call the 
calibrated GRBs HD, consisting of a sample of 109 objects, shown in Fig. HI Their redshift 
distribution covers a broad range of from 0.033 to 8.23, thus extending far beyond that of 
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FIG. 4: Distance modulus fj,{z) for the calibrated GRBs Hubble diagram made up by fitting the 
Amati correlation. 

SNela {z <~ 1.7), and including GRB 092304, the new high-z record holder of Gamma Ray 
Bursts. Here we want to use such calibrated GRBs HD to test if our cosmological model is 
able to describe the background expansion up to redshifts 2 ~ 8. In our Bayesian approach 
to model testing, we explore the parameter space through the likelihood function 



Cgrb{p) oc exp [-XGiiB(p)/2] , 



(38) 



with 



Xgrb[P) 



E 

i=l 



GRB 



(39) 



where aoRs takes into account the intrinsic scatter inherited from the scatter of GRBs 
around the Amati correlation (see Ref. 3J] and references therein), p denotes the set of 



model parameters {B and n and h in our case), and the distance modulus fi{z) is provided by 
Eq. f l28|) . The inferred confidence intervals (at 3 a) for Hq, and h, are Hq G (0.96, 1.1), 
Qrn G (0.26,0.39), and h e (0.65,0.74). We obtain Xred = 0-97 for 109 data points. In Fig. 
I5]we show the GRBs Hubble diagram with overplotted the distance modulus predicted by 
the fiducial model. It turns out that our cosmological model is fully compatible with this 
recently compiled GRBs HD. 
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FIG. 5: The calibrated GRBs Hubble diagram with overplotted the distance modulus predicted by 
the fiducial model (solid line). The full circles correspond to the GRBS data set, while the empty 
red triangles correspond to the Union2 Snela data points. 



CONSTRAINTS FROM CHANDRA X-RAY OBSERVATIONS OF LARGE RE- 
LAXED GALAXY CLUSTERS 

The matter content of the largest clusters of galaxies is expected to provide an almost 



fair sample of the matter content of the Universe (see, for instance, Refs. 49|and 35|]). The 
ratio of baryonic-to-total mass in clusters should, therefore, closely match the ratio of the 
cosmological parameters fib/^m- The baryonic mass content of clusters is dominated by the 
X-ray emitting gas, the mass of which exceeds the mass of optically luminous material by a 
factor ~ 6, with other sources of baryonic matter being negligible. The combination of robust 
measurements of the baryonic mass fraction in clusters from X-ray observations together 
with a determination of fib from other measuremets (as for instance cosmic microwave 
background (CMB) data or big-bang nucleosynthesis calculations) and a constraint on the 
Hubble constant, can therefore be used to measure and constrain the parameters which 
characterize any cosmological model. This constraint originates from the dependence of 
the /gas measurements, which derive from the observed X-ray gas temperature and density 
profiles, on the assumed distances to the clusters, /gas oc d^'^. 

To understand the origin of the /gas oc d^^^ dependence, consider a spherical region of 
observed angular radius 6, within which the mean gas mass fraction is measured. The 
physical size, R, is related to the angle 6 as R = Od^- The X-ray luminosity emitted from 
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within this region, Lx, is related to the detected flux, Fx, as Lx = indf^Fx, where c^l is 
the luminosity distance and = + zY is the angular diameter distance. Since the 

X-ray emission is primarily due to collisional processes (bremsstrahlung and line emission) 
and is optically thin, we may also write Lx oc n^V, where n is the mean number density 
of colliding gas particles and V is the volume of the emitting region, with V = 4:TT{6dA)^ /3. 
On considering the cosmo logical distance dependences, we see that n oc di^/d\^, and that 
the observed gas mass within the measurement radius Mgas oc nV oc di^d\^. The total mass, 
Mtot, determined from the X-ray data under the assumption of hydrostatic equilibrium, 
is such that M^ot oc d^- Thus, the X-ray gas mass fraction measured within angle 9 is 
/gas = Mgi^/Mtot oc di^d'j^^. The expectation from non-radiative hydrodynamical simulations 
is that for the largest {kT >z 5keV), dynamically relaxed clusters and for measurement radii 
beyond the innermost core (r >z ^2500); /gas should be approximately constant with redshift, 
the virial radius r25oo being defined as the radius of a sphere such that the mean density 
contained within it is A = 2500 times the critical density at the halo redshift. 

It is worth noting that even if the virial radius r25oo depends on the fiducial cosmological 
model, in order to determine constraints on cosmological parameters it is not necessary to 
generate /gas datasets for every cosmology of interest and compare them with the expected 
behaviour. Indeed, it is possible to fit a single fiducial /gas dataset with a model that accounts 
for the expected apparent variation in /gas(^) as the underlying cosmology is varied. Let us 



choose the ACDM reference cosmology. Following Ref. 35|], the model fitted to the reference 
ACDM data is 



/gl™(^) 







1.5 




d^iz) 





(40) 



'^"^ l + s{z) 

where dA{z) and d^^^^{z) are the angular diameter distances to the clusters in our test and 
reference cosmological models, respectively. 

In order to construct the angular diameter distance for our cosmological model, we use 
the relation between the angular diameter distance dA and the luminosity distance dr, 

di^ = {l + zfdA, (41) 

di^ being given in Eq. f p3|) . In Eq. f l40|) A takes into account the change in angle subtended 
by r25oo as the underlying cosmology is varied, and can be evaluated as 



^"l^250oj l[if(z)dA(.)r°^ ' • ^''^ 
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FIG. 6: The variation of the X-ray gas mass fraction measured within r25oo as a function of redshift 
for our model. The dashed hues, which outhne the maximum likelihood region, correspond to the 
extreme values for the parameters Og and while the gray solid line corresponds to ctg = ab = 0. 



We take the value of the slope, r], of the fgas{r/r2^m) data in the region of r25oo (as measured 
for the reference ACDM cosmology), indicated in Ref. jsSj, that is = 0.214 ± 0.022. The 
parameter 7 in Eq. fHU]) models non-thermal pressure support in the clusters. On relying 
upon hydrodynamical simulations, we take 1.0 < 7 < 1.2. The parameter s{z) = SQ^l + agz) 
in Eq. ( HOj) models instead the contribution to the baryonic mass given by stars. The factor 
b{z) = 60(1 + is the 'depletion' or 'bias' factor and describes, in a completely empirical 
way, the ratio by which the baryon fraction measured at r25oo is depleted with respect to the 
universal mean. According to Ref. |35| we choose sq = (0.16 ± 0.05)/^?o^ -0.2 < < 0.2, 
0.65 < 60 < 1.0 and —0.1 < ab < 0.1 (which corresponds to a moderate, systematic evolution 
in b{z)). The factor K in Eq. fHU]) is a calibration constant fixed to K = 1.0 ± 0.1. 

As above, we perform a Bayesian analysis, maximizing our likelihood C = exp (— 
on a grid in the space of parameters B and n, and varying all the astrophysical parameters 
appearing in Eq. (1401) . Moreover, as in Ref. (sS], we use the standard priors with f2b/i^ = 
0.0214 ±0.0020 and h = 0.72 ±0.08. We obtain, in such a way, a sort of maximum likelihood 
region, where f{r 2500) gash^'^ can vary, as shown in Fig. El The inferred region of confidence 
(at 3 a) for Hq and are (0.85, 1.05) and (0.26,0.51), respectively. 



19 



COSMOGRAPHY 



General approach 

Over the last years the cosmographic approach to cosmology gained increasing interest 
for catching as much information as possible directly from observations, retaining the min- 
imal priors of isotropy and spatial homogeneity and leaving aside any other assumptions. 
Actually, the only ingredient taken into account a priori in this approach is the FLRW line 
element obtained from kinematical requirements 

dr'^ 



ds^ = -c^dr + a\t) 



1 — kr'^ 



+ r'dQ' 



(43) 



By using this metric, it is possible to express the luminosity distance (i^ as a power series in 
the redshift parameter z, the coefficients of the expansion being functions of the scale factor 
a{t) and its higher-order derivatives. Such an expansion leads to a distance - redshift relation 
which only relies on the assumption of the Friedmann-Lemaitre-Robertson-Walker metric, 
thus being fully model independent since it does not depend on the particular form of the 
solution of cosmic equations. For this purpose, it is convenient to introduce the following 
parameters: 

1 da , , 

H = -- 44 
a dt 

1 d^a , , , 

1 (47) 



aR^ dt^ 

(48) 



These parameters are usually referred to as the Hubble, deceleration, jerk[59|], and snap 
parameters, respectively. 

Their present day values (which, as usual, we will denote with a subscript 0) can be 
used to characterize the evolutionary status of the Universe. For example, qo < denotes an 
accelerated expansion, while a change of sign of j (in an expanding universe) signals that the 
acceleration starts increasing or decreasing. Most importantly, the parameters {qo,jo,SQ} 
can be used to evaluate different distances in the universe. This can be achieved by inverting 
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the scale factor series expansion (approximated to the fourth order in t — to) of a FLRW 
metric in terms of time, given in the following equation: 

^ = 1 + ffo(t - to) - |i^o^(t - to)^ + |i^o^(t - tof + |ffo^(t - tor ■ (49) 

Therefore, one can obtain a series expansion of a distance, D{ti,tQ), travelled by a given 
photon that was emitted at ti and detected at the current epoch to, in terms of the scale factor 
or redshift, while the coefficients of the expansion are defined through the cosmographic 
parameters. Such a distance can be related to several physical quantities, for example 
the luminosity distance, the angular diameter distance and more. These quantities can be 
constrained observationally through, for example, SNela, gravitational lenses, and possibly, 
GRB data. It is worth noticing, in this respect, that since the cosmographic approach is 
based on a Taylor expansion of the scale factor, or redshift, for data of GRB at high redshift 
(above 2; = 1), it is better to use the variable y = z/{l + z), introduced in Ref. |5l|, in such 
a way that z G (0, 00) is mapped into y G (0, 1), obtaining 

ddy) = -^^[y-li^o- 3)2/' + ^ [12 - 5go + Sql - jo] y' + ^ [60 - 7jo 

-10 - 32go + lOgoJo + QQo + 21ql - ISg^ + So] y' + 0{y')} . (50) 

Application of cosmography to our model 

In this subsection, we will relate the parameters characterizing the model introduced 
in Sec. II to the cosmographic parameters {qo,jo,So}. By expanding our approximate 
luminosity distance up to the fourth order in the y-parameter, and comparing such an 
expansion with the standard expansion to the fourth order, we get the map which relates 
the parameters B and n of our model to the cosmographic parameters go, jo, so- Actually, 
we find 

_ 6B{B + 2)n-B{B + 3) + 12n-2 



jo = 18 {{B + 2)V - 54(5 + 2){B{B + 2) + 2)n^ + 9(5 + 2){B{5B + 7) + 6)n 

1 

L8(5 + 2)% 



3(5 + 1)(5(45 + 9) + 8)n + (5 + 1)^(5 + 2)) | , (52) 
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So = [S^(2(n-3)n + l)(3(n-2)n + l)(6(n- l)n + l)+ 

+ B^(n(6n(6n(n(4n(2n - 9) + 59) - 43) + 85) - 79) + 5) + 

+ 52(n(n(24n(9n(2n(2n - 7) + 19) - 112) + 851) - 140) + 9)+ 

+ B{n{72n{n{2n{S{n - 3)n + 27) - 31) + 10) - 115) + 7)+ 

+ 2(6(n - l)n + 1) (4n^ - 6n + l) (6n(2n - 1) + 1)] 



(53) 



36(5 + 2)4n6 ■ 

In order to constrain the model, we need to constrain observationally the cosmographic 
parameters by using appropriate distance indicators. Moreover, we must take care that the 
expansion of the distance related quantities in terms of {qo,jo,So) closely follows the exact 
expressions over the range probed by the data used. Taking SNela and a fiducial ACDM 
model as a test case, one has to check that the approximated luminosity distance deviates 
from the ACDM one by less than the measurement uncertainties up to 2; ~ 1.5, to avoid 
introducing any systematic bias. Since we are interested in constraining the cosmographic 
parameters, we will expand the luminosity distance D^, up to the fifth order in z which 
indeed allows us to track the ACDM expression with an error less than 1% over the full 
redshift range. To constrain the parameters {h,qo,jo,So), we define the likelihood jC{p) as 

i-Xlnela/GRB/'^) 



(2vr)- 
1 



Snela/GRB . 

2 \CsneIa/GRB\ ' 



X 



X 



X 



X 



exp 



1 f h - h phs 
2 



exp(-x|Ao/2) 

(27r)A^B.o/2|C^^^|l/2 



exp 



na. 



n 



n-n 



obs 



exp (-X^/2) 



(54) 



where 



Snela/GRB 



XSnela/GRB 



;p) = E 



i=l 



l^obsi^i) — l^thi^i, P) 



(55) 



As observational dataset we actually use the currently available observational SNIa and 
GRB Hubble Diagrams, and we set Gaussian priors on the distance from Baryon Acoustic 
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Oscillations (BAO), and the Hubble constant h. Such priors have been included in order 
to help break the degeneracies among the parameters of the cosmographic series expansion. 
In Eq. Csneia/GRB IS the Snela/GRBs diagonal covariance matrix and {hobs^cFh) = 

(0.742,0.036). The third term takes into account the constraints on = Tg^Zct) / Dy{z) 
with rs{zd) the comoving sound horizon at the drag redshift Zd (which we fix to be rs{Z(i) = 
152.6 Mpc from WMAP7) and the volume distance is defined as in Ref. 52( |: 

, n 2' 



Dviz) 



cz 



H(z) 



l + z 



1/3 



(56) 



The values of aX z = 0.20 and z = 0.35 have been estimated in Ref. 
DR7 galaxy sample so that we define Xbao ~ 

D^C^^oC with D"^ = {d°^' 



53j using the SDSS 



0.2 



Jth Jobs Jth 



^0.2; "0.35 — ."fUs) 

and Cbao the BAO covariance matrix. The next term refers to the shift parameter 5J] 55| : 

dz' 



7^ = 



M 



(57) 



'0 E{z') 

with z^, = 1090.10 the redshift of the last scattering surface. We follow again WMAP7 by 
setting (7lobs,cr-R) = (1.725,0.019). Thus we obtain the following confidence region (at 3cr) 



for the cosmographic parameters: 



v 



go -0.76 -0.25 
jo 0.28 0.33 
So -2.15 -1.18 



confirming an accelerated stage of 



the universe. 

In addition, we can investigate the possibility to use high redshift GRBs to determine 
parameters of our cosmography. As a matter of fact, we use a whole dataset containing 
both the SNIa Union2 dataset and the calibrated GRBs HD, which we call the cosmographic 
dataset. We obtain the following confidence region (at 3cr) for the cosmographic parameters: 
^ qo -0.6 -0.2^ 
jo 0.01 0.21 



In Fig. [7]we actually show the cosmographic distance modulus together 



yso -2.14 0.86 J 
with the cosmographic dataset. 

As a final step, we can use the re-parametrization of our model in terms of Hq and 

Qm exhibited in the subsection Q to construct d'^^''^"^'^^^^^'^{y,HQ,Qm,h), and perform the 

cosmographic analysis using both Supernovae and Gamma Ray Bursts data (the so-called 

cosmographic dataset), which allow us to obtain constraints (even if not yet stringent) on 

the parameters of cosmography. Actually, it turns out that 



Ho = 0.98 



+0.04 
-0.05 ' 



0.2 ±0.05 



(58) 
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FIG. 7: Distance modulus for the best-fit values of our cosmography (solid blue line) performed 
with the cosmographic dataset (empty black diamond). The red vertical line indicates the extreme 
value of the parameter y = z/{l + z) where we have SNela data. 

and then, from Eqs. (El]), we obtain the following confidence region (at 3a) for the cosmo- 



graphic parameters: 



^ qo -0.77 -0.26 ^ 
jo 0.11 0.39 



In Fig [7] we plot the distance modulus for the 



-2.2 -0.64 J 

best-fit values of our cosmography (solid blue line) performed with the cosmographic dataset 
(empty black diamond). The red vertical line indicates the extreme value of the parameter 
y = z/{l + z) where we have SNela data. 

A MORE REALISTIC DESCRIPTION 

The model of the universe adopted so far is described by an exact solution of the dynami- 
cal equations, whose arbitrary parameters are determined by specifying the initial conditions. 
In this section we now consider more realistically the inclusion of radiation, too, into such 
a model. In this case the dynamical equations, as far as we know, do not have analytical 
solutions, and therefore we will rely on numerical computations. 

The dynamical field equations become 

- 1 - i fl + '-^^^-^ - 1)"'--"'-' - ° • (^^) 

,G_5,|,3,^G.|(-6^;-2A + 2G§)^0. (60) 
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We have also to consider the Hamiltonian constraint [2 

A fiG^ SttG 



a 



2 



3 6G2 3 



(D„a-3^-+D,,,a-3>-)=0. (61) 



It turns out that these equations assume a simpler form when, instead of t as an independent 
variable, one uses a{t) - the scale factor. We introduce a new independent variable by 
u = log(l + z) = — log(^), where ao is the present value of the scale factor (fixed at 
Oq = 1) and z is the redshift. The equations can now be written in the form 

, / 2 H'(u)\ A(u) 2 ,/G'(u)\^ Stt^, , , 

G(.) = G(.)^3 + ^ + -^J+^^-3 + ^^-^^^^^J , (63) 

8^P^^^a-37..d + D^a-3^™) + 4 , , 

^'(«) = ' t;;-^ , 64 

G(u) 

where ' = ^, 7r.ad = |, and 7^ = 1. 

Moreover, we have to remember that above, in Eq. (!63|) . we explicitly used the dependence 
of A on G, as given by applying the Noether Symmetry Approach to the case of a matter- 



dominated universe (without any radiation content) 60]: A(u) = WG^-^'^ (u). Let us note 



that Eqs. ( I62l) and ( l64l) can be rewritten in a Friedmann-like form 

, / 2H'(u)\ AttG , , 
^'(«) (^1 - = (PA,G +p™d) , (65) 

tt2/ \ SnGiu) , . 

(^) = ^ (PA,G + Prad + Pm) , (66) 

where we define 

. - ^ I Hju)^ (3 - t) 

^'^'^ = 8^ + ' (^^^ 

_ H^3-f)uHGf A 
^^'^ = 4^ ' (6^^ 

and then construct w = Substituting Eq. f l64|) in Eq. f l63|) . we can numerically 

Pa,g 

solve the system of differential equations characterizing cosmology, by specifying the initial 
condition at m = 70, for example, and assuming that G{70), and G'{70) have the same values 
as in the case without radiation. As can be seen in Fig. |H1 the presence of radiation affects 
the evolution of the Q parameters. As expected, at the initial time, again when u = 30, 
radiation dominates the expansion rate of the universe, with the A-term and matter being 
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FIG. 8: parameters as functions of u in the universe filled in with matter, radiation and the 
A-term. is marked in blue, Vl^ in red and f^m in green. 




FIG. 9: The effective equation of state u; as a function of u. 

subdominant, at a redshift z of about 5000; the energy densities of matter and radiation 
become equal, and, for a relatively short period, the universe becomes matter dominated 
until, at a redshift of about 1, the A-term starts dominating the expansion rate of the 
universe. In Fig. |9]we plot the behaviour of w as a function of u: we can therefore see that 
it behaves like stiff matter {w = |) up to u ^ 10, when there is the beginning of a transition 
towards a superquintessence behavior with w < —1. 

It is worth noting that, since the initial conditions needed to numerically integrate the 
Eqs. fl65|) . fl66|) and f l67|) . have been set with the only constraint of obtaining a well-behaved 
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evolution for the Q parameters, the values w{0), ftA{0), fir(O), and $7^(0), cannot be directly 
compared with the ones indicated in the literature (see for instance Ref. |33|), because 
they do not result from a fitting procedure on observational datasets. Moreover, as far 
as the equation of state of dark energy is concerned, the best-fit values also depend on the 
mathematical law assumed for it. For instance, if we adopt models of dark energy as potential 
energy of some, as yet undiscovered, scalar field, we cannot obtain superquintessential values 
w < —1, unless we consider phantom fields or non-minimal coupling terms. On the other 
hand, such values are fully acceptable when we look at the confidence regions obtained by 
fitting the data with a parametrized form of the equation of state; however, this is not the 
same as saying that the scalar field is ruled out by the statistics. 



CONCLUSIONS 



We have shown that one can build a matter-dominated cosmological model with variable 
Newton parameter and variable cosmological term which is compatible with the more up- 
dated observations of type la supernovae, the gamma ray bursts Hubble diagram, and the 
gas mass fraction in X-ray luminous galaxy clusters. Moreover, we have applied to such a 
cosmological model a cosmographic approach, which can help in selecting realistic models 
without a priori choices that can be questionable. In performing our cosmographic analysis 
we set Gaussian priors on the distance from Baryon Acoustic Oscillations (BAO), and the 
Hubble constant h. Such priors have been included in order to help break the degeneracies 
among the parameters of the cosmographic series expansion. A more realistic approach, 
considering the inclusion of a radiation component, seems also possible, even if it has to be 
worked out only numerically. 

Some questions remain, however, unsolved; for example, it is not clear enough what can 
be definitely said about the effects of G, which might point out to the violation of the strong 
equivalence principle, but not necessarily of the Einstein equivalence principle. Here we have 
simply discussed the observational bounds on the possible variations of the gravitational 
constant in the early universe, considering the best limit (at 3a), 5G = O.OQlojg, obtained 
in the literature. It turned out that our model can satisfy such a best fit limit, provided 
that the n parameter is appropriately chosen (the role of B is only marginal with respect to 
this strong test) to be compatible with the basic cosmological observations. We postpone 
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to a forthcoming paper the detailed analysis of the the dependence of the various elemental 
abundances on the time variation of G for our model. We expect that the BBN constraints 
are the main tool to restrict the allowed region of the model parameters, so representing a 
necessary step in the feasibility study of the model itself. 

The concrete possibility to generate acceptably a previous radiation-dominated regime in 
the framework of Renormalization-Group inspired cosmology has still to be proved. Unfor- 
tunately, the procedure of the Noether Symmetry Approach does not work with a Lagrangian 
where the matter term is Lm = -0(1 — a^^) (being 7 = 1 for dust and 7 = 4/3 for radiation). 
However, as we have seen, a numerical integration of the equations derived by using such an 
Lmi implies that, after inflation, such equations do give rise to a period of radiation domi- 
nance followed by one of matter dominance. Later, a period of acceleration only depending 
on variability of A becomes the dominant one. 

As far as we can see, we think that the work here presented may be considered as 
another important step towards a more accurate confrontation of variable-G cosmologies 
with modern observations. Even if, as said in the Introduction, such a comparison has begun 
earlier with many of the papers cited in the references, they all indeed only contribute to 
show a first feasibility of some of such cosmological models, not only from a theoretical 
but also from an observational point of view. Much work still remains to be done but the 
years to come will hopefully help in further restricting the family of these and other possibly 
viable cosmological models, mainly from the observational point of view. Actually, in a 
forthcoming paper we are going to investigate the behavior of density perturbations during 
the matter dominated stage. 

Acknowledgments 

G. Esposito and C. Rubano are grateful to the Dipartimento di Scienze Fisiche of Federico 
II University, Naples, for hospitality and support. 



* E-mail at: lester@na.infn.itl 
[1] Davis, T. M., et al: Astroph. J. 666, 716 (2007) 
[2] Ostriker, J. P.: Ann. Rev. Astron. Astrophys. 31, 689 (1993) 



28 



[3] Lusanna, L.: "Dark matter as a relativistic inertial effect in Einstein canonical gravity?" 
(2010) (|arXiv:1011.2 908[gr-qc]) 

[4] Durrer,R., Maartens, R.: Gen. Relat. Grav. 40, 301 (2008) 

[5] Capozziello, S., Francaviglia, M.: Gen. Relat. Grav., 40, 357 (2008) 

[6] Capozziello, S., Faraoni, V.: Beyond Einstein Gravity, Springer, Dordrecht, 2011 

[7] Renter, M., Wetterich, C.: Nucl. Phys. B 427, 291 (1994) 

[8] Renter, M.: Phys. Rev. D 57, 971 (1998) 

[9] Souma, W.: Prog. Theor. Phys. 102, 181 (1999) 
[10] Wetterich, C.: Int. J. Mod. Phys. A 16, 1951 (2001) 
[11] Berges, J., Tetradis, N., Wetterich, C.: Phys. Rep. 363, 223 (2002) 
[12] Lauscher, O., Renter, M.: Phys. Rev. D 65, 025013 (2002) 
[13] Renter, M., Saueressig, F.: Phys. Rev. D 65, 065016 (2002) 
[14] Lauscher, O., Renter, M.: Class. Quantum Grav. 19, 483 (2002) 
[15] Babic, A., Guberina, B., Horvat, R., Stefancic, H.: Phys. Rev. D 71, 124041 (2005) 
[16] Codello, A., Percacci, R.: Phys. Rev. D 97, 221301 (2006) 
[17] Bonanno, A., Renter, M.: JHEP 02, 035 (2005) 
[18] Niedermaier, M.: Nucl. Phys. B 673, 131 (2003) 
[19] Niedermaier, M.: JHEP 12, 066 (2002) 

[20] Forgacs, P., Niedermaier, M.: "A fixed point for truncated quantum Einstein gravity" (2002) 

( |hep-th/ 0207028j) 
[21] Bonanno, A., Reuter, M.: Phys. Rev. D 65, 043508 (2002) 
[22] Bonanno, A., Reuter, M.: Phys. Rev. D 60, 084011 (1999) 
[23] Bonanno, A., Reuter, M.: Phys. Rev. D 62, 043008 (2000) 
[24] Tsamis, N. C, Woodard R. P.: Phys. Lett. B 301, 351 (2003) 
[25] Bonanno, A., Reuter, M.: Phys. Lett. B 527, 9 (2002) 
[26] Bentivegna, E., Bonanno, A., Reuter, M.: JCAP 01, 001 (2004) 
[27] Bonanno A., Esposito, G., Rubano, C: Class. Quantum Grav. 21, 5005 (2004) 
[28] de Ritis, R., Marmo, G., Platania, G., Rubano, C, Scudellaro, P., Stornaiolo, C: Phys. Rev. 

D 42, 1091 (1990) 

[29] Capozziello, S., de Ritis, R., Rubano, C, Scudellaro, P.: Rivista del Nuovo Cimento 19(4), 1 
(1996) 

29 



[30] Bonanno, A., Esposito, G., Rubano, C, Scudellaro, P.: Class. Quantum Grav. 24, 1443 (2007) 
[31] Bonanno, A., Esposito, G., Rubano, C., Scudellaro, P.: Int. J. Geom. Meth. Mod. Phys. 5, 

329 (2008) 

[32] Bonanno, A., Esposito, G., Rubano, C., Scudellaro, P.: Gen. Relat. Grav. 39, 189 (2007) 
[33] Amanullah, R., et al.: Astroph. J. 716, 712 (2010) 

[34] Demianski, M., Piedipalumbo, E., Rubano, C.: Mon. Not. R. Astron. Soc. 411, 1213 (2011) 
[35] Allen, S. W., Rapetti, D. A., Schmidt, R. W., Ebeling, H., Morris, G., Fabian, A.C.: Mon. 

Not. R. Astron. Soc. 383, 879 (2008) 
[36] Capozziello, S., dc Ritis, R.: Gen. Relat. Grav. 29, 1425 (1997) 
[37] Capozziello, S., de Ritis, R., Marino, A. A.: Gen. Relat. Grav. 30, 1247 (1998) 
[38] Bonanno, A., Esposito, G., Rubano, C, Scudellaro, P.: Class. Quantum Grav. 23, 3103 (2006) 
[39] Bambi, C, Giannotti, M., Villante, F. L.: Phys. Rev. D 71, 123524 (2005) 
[40] Pitjeva, E. V.: Astron. Lett. 31, 340 (2005) 

[41] Jarosik, N., Bennett, C. L., Dunkley, J., Gold, B., Greason, M. R.: Astroph. J. Suppl. S. 192, 
14 (2011) 

[42] Schmidt, B. P., et al.: Astroph. J. 507, 46 (1998) 
[43] Perlmutter, S., et al.: Astroph. J. 517, 565 (1999) 

[44] Gaztanaga, E., Garcia-Berro, E., Isern, J., Bravo, E., Dominguez, I.: Phys. Rev. D 65, 023506 
(2002) 

[45] Uzan, J. Ph.: Rev. Mod. Phys 75, 403 (2003) 

[46] Rubano, C, Scudellaro, P., Piedipalumbo, E., Capozziello, S., Capone, M.: Phys. Rev. D 69, 

103510 (2004) 
[47] Schaefer, B. E.: Astroph. J. 660, 16 (2007) 

[48] Demianski, M., Piedipalumbo, E.: "Standardizing the GRBs with the Ep^i - Eis^ relation: the 
updated Hubble diagram and implications for cosmography" (2011), submitted to Mon. Not. 
R. Astron. Soc. 

[49] White, S. D. M., Navarro, J. F., Evrard, A. E., Frenk, C. S.: Nature 366, 429 (1993) 

[50] Sahni, V., Saini, T. D., Starobinsky, A. A., Alam, U.: JETPL 77, 201 (2003) 

[51] Chevallier, M., Polarski, D.: Int. J. Mod. Phys. D 10, 213 (2001) 

[52] Eisenstein, D.J., et al.: Astroph. J. 633, 560 (2005) 

[53] Percival, et al.: Mon. Not. R. Astron. Soc. 401, 2148 (2010) 

30 



[54] Bond, J.R., Efstathiou, G., Tegmark, M.: Mon. Not. R. Astron. Soc. 291, L33 (1997) 
[55] Efstathiou, G., Bond, J.R.: Mon. Not. R. Astron. Soc. 304, 75 (1999) 

[56] It is worth noting that such a constraint on the parameters pi of our model is not strictly 



311 ] ■ as an instrument to fix somehow the 



necessary, but it has been used, already in Ref. 
reachable regions in the space of parameters (at least for what it concerns n). This means that 
one could perform our analysis leaving the pi parameters, and then n, completely free 

[57] This limit refers to the value of G when the temperature of the universe is 0.02 < T < 0.2 
MeV (i.e. during and immediately after the d-bottleneck epoch) and is consistent with the 
standard assumption that G has lightly varied during the evolution of the universe. 

[58] Let us note that in principle it could be better to use the re-parametrization described above 
(i.e. the dependence of B and n on Hq and Jim provided by Eqs. (j23p and (j24p ) and to 
maximize the likelihood with respect to the parameters h, Hq and since they have a 
brighter physical meaning and their space of parameters is easily reconstructed as discussed 
above. However, the mathematical form of B{HQ,0,rn) and n{HQ,Qrn), when inserted in the 
luminosity distance and in the correction term AMq, is such that they make the numerical 
calculations too hard and time consuming. Thus, we prefer to work in the space of parameters 
h, B and n. Anyway, we will work in the space of parameters h, Hq and later, when we 
implement the cosmographic approach. 

between different models was also proposed in 



5C] 



[59] The use of the jerk parameter to discriminate 

the context of the statefinder parametrization 
[60] We actually remember that in the presence of radiation it is not possible to solve analytically 

the equations by defining the existence of the Noether symmetry vector field. 



31 



